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^ ■ Three ways are compared of dealing with rescattering effects in B"^ 

Q"^ ■ K^TT^, in order to determine the weak phase 7 from these processes and 

r-| ■ -B^ — ^ K^7r°. We find that neglecting these contributions altogether may 

Qh! involve sizable errors in 7, depending on the rescattering amplitude and 

Q^l on the value of a certain measurable strong phase. We show that an 

attempt to eliminate these effects by using the charge-averaged rate of 
B^ -^ K^K^ suffers from a large theoretical error due to SU(3) breaking, 
which may be resolved when using also the processes B^ -^ t^^Vs- 



PACS codes: 12.15.Hh, 12.15.Ji, 13.25.Hw, 14.40.Nd 



1 Introduction 

The weak phase 7 = Arg(V^*^) is presently the least well known quantity among 
the four parameters (three angles and a phase) of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. Its determination, which is regarded to be more difficult than that of 
the other two angles of the CKM unitarity triangle ^, 0, can provide a crucial test of 
the CKM mechanism for CP violation in the Standard Model. Several methods have 
been proposed to determine 7 from hadronic two-body B decays. The methods which 
seem to be experimentally most feasible in the near future are based on applications 
of SU(3) flavor symmetry in B decays into two light charmless pseudoscalars [Q]. 
These methods involve certain theoretical uncertainties, which are expected to be 
reduced when more data become available and when better theoretical understanding 
of hadronic B decays is achieved. 

In a first paper in a series, Gronau, London and Rosner (GLR) Q proposed to 
extract 7 by combining decay rate measurements of B~^ -^ Ktt, B^ -^ vrvr with their 
charge-conjugates. SU(3) breaking, occuring in a relation between B ^ titi I = 2 
and B — > Kn I = 3/2 amplitudes, was introduced through a factor fx/ fw when 
assuming that these amplitudes factorize. In its original version, suggested before 
the observation of the heavy top quark, the method of Ref. [Q neglected electroweak 
penguin (EWP) contributions and certain rescattering effects. Subsequently, model- 
calculations showed that due to the heavy top quark the neglected EWP terms were 
significant 0; and recently these terms were related by SU(3) to the B —>■ Kn I = 3/2 
current-current amplitudes P, 0] . This led to a modification [H of the GLR method, 
to be referred to as the GLRN method, which in the limit of fiavor SU(3) symmetry 
includes EWP effects in a model-independent way. Corrections from SU(3) breaking, 
affecting the relation between EWP terms and current-current terms, were argued to 
be small ||, |. 

Assuming that the above SU(3) breaking effects are indeed under control, there is 
still an uncertainty due to rescattering effects. To determine 7 from the above rates, 
one takes the B'^ — * K^tt^ amplitude to be pure penguin, involving no term with 
weak phase 7. This assumption, which neglects quark annihilation and rescattering 
contributions from charmless intermediate states, was challenged by a large number 
of authors |T^. Several authors proposed ways of controlling rescattering effects 
in B^ -^ K^TT^ by relating them through SU(3) to the much enhanced effects in 
B^ -^ K^K^ [0, ig 0] (see also 0, |TB|, |T§). The charge-averaged rate of the 



latter processes can be used to set an upper limit on the rescattering amphtude in 
B^ -^ K^TT^. While present limits are at the level of 20 — 30% of the dominant 
penguin amplitude [^ Q (depending somewhat on the value of 7) , they are expected 
to be improved in the future. The smaller the rescattering amplitude is, the more 
precisely can 7 be determined from the GLRN method. A recent demonstration [0], 



based on a few possible rate measurements, seems to show that if the rescattering 
amphtude is an order of magnitude smaller than the dominant penguin amplitude in 
B^ -^ K^TT^, the uncertainty in 7 is only about 5 degrees. 

In the present Letter we reexamine in detail the uncertainty in 7 due to rescat- 
tering effects. Using a geometrical interpretation for the extraction of 7, we perform 
in Section 2 numerical simulations which cover the entire parameter space of the two 
relevant strong phases, the rescattering phase (pA and the relative phase between 
I = 3/2 current-current and penguin amplitudes. We find that, contrary to the 
demonstration made in P, a 10% rescattering amplitude leads to an uncertainty in 
7 as large as about 14° around ~ 90°. For certain singular cases no solution can be 
found for 7. We show that can be determined rather precisely from the B^ — >■ Kir 
rate measurements |^, which could reduce substantially the error in 7 if values far 
apart from = 90° were found. 



It has been suggested |T3| to go one step beyond setting limits on rescattering 
contributions in A{B^ — >■ K^n^) and to completely eliminate them by using the 
charge- averaged rate measurement of B^ — > K^K^. Applying our geometrical for- 
mulation, we will show in Section 3 that the resulting determination of 7 is unstable 
under SU(3) breaking which can introduce very large uncertainties in 7 . 

Finally, in order to overcome these uncertainties, we have recently proposed to 
use in addition to B^ -^ K^K^ also the processes B^ — >■ vr^r^g Jl^- Although this 



may be considered an academic exercise, mainly due to complicating rj — rj' mixing 
effects, we will examine in Section 4 the precision of this method. We will show that, 
when neglecting 77 — 77' mixing, the theoretical error in 7 is reduced to a few degrees. 
We conclude in Section 5. An algebraic condition, used in Section 3 to eliminate 
rescattering effects by B^ -^ K^K^ decays, is derived in an Appendix. 

2 Rescattering uncertainty in 7 from B^ -^ Ktt 

The amplitudes for charged B decays can be parameterized in terms of graphical 
contributions representing SU(3) amplitudes (we use the notations of [0]): 

A{B+^K%+) = \X'^Y'^{A + Puc) + xl'\Pct + Pi'^), (1) 

V2A{B+ ^ K+tt') = |A(^)|e*^(-T - C - A - P^,) + Xi'\-P,t + V2PD , (2) 

\/2A(5+ -> 7r+7r°) = \Xi'^\e'^-T-C) , (3) 

where X^ = V^iiXqiq are the corresponding CKM factors. These amplitudes satisfy a 
triangle relation P, p| 



-J7 



V2A(5+ ^ i^+7r°) + A(5+ ^ ir°7r+) = V2f„|A(fi+ ^ 7r+7r°)|e*(^+«) (l - bEwe 

(4) 



Here we denote f„ = (/i^//^)A/(l - AV2) ~ 0.28, 5ew = -(3/2)|A^'V^l'^l«^ - 0.66 
(k = (c9+Cio)/(ci+C2) = —8.8-10"^), while ^ is an unknown strong phase. The second 
term in the brackets represents the sum of EWP contributions to the amphtudes on 
the left-hand side [^ |^. The factor fx/ fn accounts for factorizable SU(3) breaking 
effects. 

The relation (^), together with its charge-conjugate counterpart, written for A{B -^ 
f) = e^"^"* A{B -^ /), are represented graphically by the two triangles OAA' and OBB' 
in Fig. 1. Here all amplitudes are divided by a common factor A = \/2fu\A{B^ -^ 
7r"'"7r'')|e**-'^"'"^\ such that the horizontal line 01 is of unit length and the radius of the 
circle is 5ew- Four of the sides of the two triangles are given by 

1 |A(g+^ir°7r+)| _ 1 |A(E+^ir+7r°)| 

1 \A{B- ^ K'^'K-)\ __ _ l\A{B- -^K~t:^)\ 



v/2fJA(fi+^7r+7rO)| ' -- f, |A(5+^7r+7rO)| ■ 

The relative orientation of the two triangles depends on 7 and is not determined from 
measurements of the sides alone. Assuming that the rescattering amplitude with weak 
phase 7 in B^ — > K^ir^ can he neglected, one takes the amplitude (|l]) to be given 
approximately by the second (penguin) term [^, ^], which implies OB = e^^^OA in 
Fig. 1. In this approximation, the weak phase 7 is determined by requiring that the 
angle (27) between OA and OB is equal to the angle (27) at the center of the circle 

i- 

In order to study the precision of determining in this way the phase 7 as function 

of the rescattering contribution which is being neglected, let us rewrite (0) in the form 



A{B+ ^ K\+) = -V,Al--\ p(l + eAC^^^e'^) , p^Pct + P^"^ (6) 




where e^ measures the magnitude of rescattering effects. In Fig. 1 the magnitude of 
these effects has a simple geometrical interpretation in terms of the distance of the 
point Y from the origin O, e^ = ll^OI/iy^l, where YO and Y A are the two com- 
ponents in the B^ —>■ K^ii'^ amplitude carrying weak phases 7 and zero, respectively 

YO = \Xi'^\e''[{A + P„,) - p]/A , YA = V,, (l - ^) p/A . (7) 

The rescattering phase 0a is given by 0a =ATg{YO/YZ), where Z is any point on the 
line bisecting the angle AYB. A second strong phase which affects the determination 
of 7 is 0, the relative strong phase between the penguin amphtude p and the I = 3/2 
current-current amplitude T + C. In Fig. 1 this phase is given by =ATg{YZ/OI). 



Let us now investigate the dependence of the error in 7 when neglecting rescat- 
tering on the relevant hadronic parameters. Our procedure will be as follows. First 
we generate a set of amplitudes based on the geometry of Fig. 1 and on given values 
of the parameters 7, e, e^i, (J)a and 0; then we solve the equation cos 27 = cos{BOA) 
and compare the output value of 7 with its input value. Here e is given in terms of 
the ratio of charge-averaged branching ratios [^ § 



A / 



K 



1 - AV2 /, \ B{B^ -^ K^ 
The geometrical construction in Fig. 1 is described by 



25(B± ^ vrivrO) ^^^ 



7r±l 



g«(</'-7) 



YA = 01 , OY = eAC^^'^^+^^A , (9) 

eJl + 2eA cos 0^ cos "y + e\ 



implying a rate asymmetry between B^ -^ K^tt^ and B -^ K^n . 

For illustration, we take 7 = 76°, e = 0.24 0], e^ = 0.1 (which is a reasonable 



guess JTT], |T2[), and we vary and 0a in the range 0° < < 180°, -90° < 0a < 270°. 
The results of a search for solutions in the interval 65° < 7 < 90° are presented in 
Fig. 2 which displays a twofold ambiguity. Fig. 2(a) shows the solution as function 
of 0A for two values of 0, = 60° and = 90°. Whereas for 0a = 90° the solution 
is very close to the input value, the deviation becomes maximal for 0a = 0°, 180°. 
This agrees with the geometry of Fig. 1, in which the largest rescattering effects are 
expected when YO is parallel or anti-parallel to the line bisecting the angle BY A. 

In a second plot. Fig. 2(b), we fix 0a = 0° and vary over its entire range, which 
illustrates the maximal rescattering effect. We find two branches of the solution for 7, 
both of which deviate strongly from the input value 7 = 76° for values of around 90°. 
At = 90° there is no solution for eA = 0.1 in the considered interval. We checked 
that the solution is restored and approaches the input value as the magnitude of eA 
decreases to zero, as it should. Thus, the uncertainty in 7, seen both in Fig. 2(a) and 
Fig. 2(b) at 0a = 0° and around = 90°, is about 14°. It can even be worse in the 
singular cases where no solution for 7 can be found. 

A variant of this method for determining 7, proposed recently in 0, was formu- 
lated in terms of two quantities R^ and A defined by 

_ B{B^ -. K^n^) 
* ~ 2B{B^ ^ ir±7rO) ' ^ ' 

~ _ B{B+ -^ K+TT°) - B{B- -^ K-7T°) B{B+ -^ K°Tr+) - B{B- -^ K^tt-) 



B{B^ -^ KH^) 2B{B^ -^ KH^) 

These quantities do not contain C(eA) terms; their dependence on the rescattering 
parameter ea appears only at order 0{eeA)- Therefore, it was argued in 0, the 



determination of 7, by setting e^i = in the expressions for /2* and A, is insensitive 
to rescattering effects. This procedure gives two equations for 7 and which can be 
solved simultaneously from i?^, and A. Using two pairs of input values for {R^:,A) 
(corresponding to a restricted range for (pA and (j)) seemed to indicate that the error 
in 7 for ca = 0.08 is only about 5°. (the relations between the parameters used in [Q 
and ours are = —0, ?] = 0^ + tt, €3/2 = e and e^ = eyi)- 

In Fig. 3 we show the results of such an analysis carried out for the entire parameter 
space of (pA and 0. Whereas the angle can be recovered with small errors, the results 
for 7 show the same large rescattering effects for values of around 90° as in Fig. 2. 
(A slight improvement is the absence of a discrete ambiguity in the value of 7.) 
These results show that the large deviation of 7 from its physical value for = 90° 
is a general phenomenon, common to all variants of this methods. Some information 
about the size of the expected error can be obtained by first determining 0. Values 
not too close to 90° would be an indication for a small error. 

3 Eliminating rescattering by B"^ -^ K^K^ 

The amplitude for B^ —>■ K^K^ is obtained from A{B^ -^ K^tt^) in (|l]) by a f/-spin 
rotation ITll 



AiB+^K+K') = \Xlf^\e'^{A + P^,) + \X[''^\e~'('{P,t + Pi'^) . (11) 



In the limit of SU(3) symmetry the amplitudes in (|TT]) are exactly the same as those 
appearing in (|I]). In Fig. 4 A{B^ -^ K^K^), scaled by the factor A/(l - \^ /2) (and 
divided by A as in Fig. 1), is given by the line OC and its charge-conjugate is given 
by OD. We have shown in |]l^ that knowledge of these two amplitudes allows one to 
completely eliminate the rescattering contribution A + Puc from the determination of 
7. This is achieved by effectively replacing in the GLRN method the origin O by the 
intersection Y of the lines AC and BD. 7 is determined by requiring that the angle 
(27) between YA and YB is equal to the angle (27) at the center of the circle. 

The amplitude (|TT|) can be decomposed into two terms carrying definite weak 
phases in form very similar to (|^), 

^ A{B^ ^ K^K') = -V^ {l-^\p f-7^^^ + eAc'^^eA , (12) 



1 - A2/2 



The ratio \CY\/\AY\ = AV(1 - X^/2f implies that the triangle AYB is about 25 
times larger than the triangle CYD. This will result in a large uncertainty in 7 also 
when the equality between the corresponding terms in B^ -^ K^n^ and B^ -^ K~^K^ 
amplitudes involves relatively small SU(3) violation. 



The geometrical construction by which rescattering amphtudes can be completely 
eliminated in the SU(3) limit consists of three steps. (See Fig. 4. For an alternative 
suggestion, see [jl3[-) 

a) Determine the position of the point F as a function of the variable angle 27 
and the decay rates of B^ -^ Kir and B~^ — >■ tt+tt'^. The point Y is chosen on the 
mid-perpendicular of AB such that the equality of the angles marked 27 is preserved 
for any value of 7. 

b) Draw two circles of radii A/(l — A^/2)|yl(i?^ — > K^K^)\ centered at the origin 
O (dashed-dotted circles in Fig. 4). The intersections of the lines AY and BY with 
these circles determine C and D respectively (up to a two- fold ambiguity), again as 
functions of 7. 

c) The physical value of 7 is determined by the requirement \AC\ = \BD\ [|1^]. 
This condition on 7 can be formulated in an algebraic form, showing that only the 
charge-averaged rate of B^ -^ K^K'^ is needed. The condition is given by Eq. (^5|) 
in the Appendix. 

Let us examine the precision of this method for e^ = 0.1 at ~ 90°, for which the 
simpler method of Sec. 2 receives large rescattering corrections. In Fig. 5(a) we show 
the left-hand side of Eq. (|r^) as a function of variable 7 at = 90° for several values 
of (pA- The value of 7 is obtained from the condition that the left-hand side of this 
equation vanishes. In the absence of SU(3) breaking this method reproduces precisely 
the physical value of 7 (7 = 76°) for all values of 0^- However SU(3) breaking effects 
can become important, to the point of completely spoiling this method. We simulate 
these effects by taking the amplitudes p and a = A+Puc—p in B^ -^ K^K^ (Eq. (pT])) 
to differ by at most 30% from those in B^ -^ K^ti"^ (Eq. ([l|)). This expands the lines 
of Fig. 5(a) into bands of finite width, which give a range for the output value of 7. 

In Fig. 5(b) we show the effects of SU(3) breaking on the determination of 7 as 
function of 0a for (f) = 90°. We see that for values of |0a| larger than about 25° 
the error on 7 is quite large. Thus, we conclude that for certain values of the strong 
phases the determination of 7 using this method is unstable under SU(3) breaking in 
the relation between B+ -^ K^7r+ and B+ -^ K+K^. 

4 The use of B^ -^ tt^Vs 



In Ref. ||T3 we proposed to use in addition to B^ -^ K^K^ also B^ — >■ 7r^?78 and 



their charge-conjugates. Writing 

A{B+ ^ TT+T^s) = |A(f)|e^^(-r - C - 2A - 2P„,) + |Ai')|e-^^(-P,, + P^^^) , (13) 



we find the triangle relation 

A{B+ ^ K^K') + \/| ^(5+ - vr+r^s) = ^^(5+ - vr+7r°) . (14) 

This relation and its charge-conjugate provide another condition which determines 
the positions of the points C and D. As in Section 3, the phase 7 is determined by the 
equation cos{BY A) = cos 27, where the point Y is fixed by the intersection of the lines 
AC and BD. General considerations, based on the relative sizes of the amplitudes 
involved, suggest that this method is relatively insensitive to SU(3) breaking effects 



17 



We illustrate this in Fig. 6 where we show on the same plot the two sides of the 
equation cos{BYA) = cos 27 as functions of the variable 7. As in method of Section 3, 
SU(3) breaking is simulated by taking the penguin (p) and annihilation (a) amplitudes 
in B^ —>■ K^K^ to differ by at most 30% (separately for their real and imaginary 
parts) from those in B^ -^ K^t^"^. The latter are used to construct the positions of 
the points C and D. In the example of Fig. 6 we take e^ = 0.1, = 90°, (pA = 45°, for 
which the two methods described in Sections 2 and 3 were shown to lead to large errors 
in 7. For an input value 7 = 76°, the output is given by the range 74° < 7 < 78°, 
obtained by the intersection of the solid line with the band formed by the diamond 
points. We see that the error in 7 due to SU(3) breaking is less than ±2°, which 



confirms the general arguments of [0. This scheme, or rather its analogous version 
using B'^ and Bs decay Jl^ , may prove useful for a determination of 7 in case that 
the strong phases (0, 0^) turn out to have values which preclude the use of the two 
simpler methods. 

5 Conclusion 

We compared three ways of dealing with rescattering effects in 5^ -^ K'^tt^, in 
order to achieve a precise determination for the weak phase 7 from these processes 
and B^ — >■ K^tt^. In the simplest GLRN method we find that large errors in 7 
are possible for a particular range of the strong phases, ~ 90°, even when the 
rescattering term is only at a level of 10%. B^ and B~ decay rate measurements 
are expected to provide rather precise information on 0. Small errors in 7 would 
be implied if turns out to be far away from 90°. The second method, in which 
rescattering effects can be completely eliminated in the SU(3) limit by using also the 
charge-averaged B^ —>■ K^K^ rate, suffers from a sizable uncertainty due to SU(3) 
breaking. These uncertainties would be resolved in an ideal world, where B^ —>■ vr^r/g 
can be measured, or alternatively by using corresponding B^ and Bs decays. 



6 Appendix 

The weak angle 7 is fixed in the method described in Sec. 3 by the condition \AC\ = 
\BD\, or equivalently \YC\ = \YD\. Exphcitly, this can be written after some algebra 
as an equation in 7 

2(1 - xofY'' + 2a;o(l - xo)Y ■(1+3) + xl{xl^ + xl_) - (y^ + y\) = , (15) 

where xq is defined as the ratio of two CP rate differences 

xo = {yl, - y\)l{xl^ - il.) s^) -AV(1 - \^/2f . (16) 

Here 

1 MAiB+^K+K^)\ __ 1 U \A{B--^K-K^)\ 

^+° ^JK \A{B+^iT+iT^)\ ' ^-° ^/^ |A(i?+^7r+7rO)| ^ ' 

obey an SU(3) relation with the amplitudes (|) of B"^ — > Kvr^ 

A^ 

This implies that CP rate differences in B^ -^ K^ti^ and B^ -^ K^K^ are equal 



and of opposite sign [13|. We see that in the SU(3) limit the condition (15), which 
eliminates rescattering effects, requires only a measurement of the charge-averaged 
rate of B^ — > K'^K^ and not the CP asymmetry in these processes |]13 . 



To prove (|T^), let us consider two lines AY and BY cutting two circles of radii -Ri, 
i?2 (centered at the origin) at points C and D respectively. The intersection points 
can be written as C = Y + Xi{A — Y) and D = Y + X2{B — Y), where Xi,X2 are 
solutions of the equations 

{A - Yfxl + 2xiY ■{A-Y) + {Y^ - Rl) = ^ (19) 

{B - Yfxl + 2x2!^ ■ (5 - f) + (f^ - i?2) = . (20) 

The condition \YC\ = \YD\ is equivalent to requiring that these two equations have 
a common solution xi = X2- Obviously, if such a solution exists, it is given by 

xo = ^^^-^^ = 5^^, (21) 

2Y-{A-B) A^-B^ 

where we used the equality {A — Yf = {B — Yf. Taking the sum of (0) and ( pOD 
with the value (0) for x leads immediately to the condition (|15|). 
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Figure 1: Relative orientation of B^ -^ Kn amplitude triangles. 
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(b) 



Figure 2: The weak phase 7 is obtained as the solution to the equation cos(27) = 
cos{BOA). (a) - the dependence of the solution on 0^, for two values of = 60° 
and = 90°; (b) - the dependence of the solution on 0, for 0^ = 0°. (both graphs 
correspond to e^ = 0.1,7 = 76°) 




80 100 
phi 



(b) 



Figure 3: (a) - the weak phase 7 extracted from the method using the parameters 
(/?*, A), as a function of the strong phase for several values of 0^ {^a = 0.1). The 
horizontal line shows the assumed physical value of 7 = 76°. (b) - the strong phase 
can be reconstructed using the [R^,, A) data. 
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Se- 



,w 



Figure 4: Geometric construction for the method described in Sec. 3. C and D denote 
the intersection points of the hues AY and BY determined as explained in the text, 
with the two circles of radii given by 1^4(5^ -^ K^K^)\. 



phi=90 




100 150 200 
phiA 



(b) 



Figure 5: (a) - The left-hand of Eq. (|T3p as a function of variable 7 for </> = 90° and 
for different values of (pA- AH these curves intersect at 7 = 76°, which is the assumed 
physical value, (b) - SU(3) breaking effects introduce an error on the extracted value 
of 7, here shown as function of 0^ at = 90°. 
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phi=90, phiA=45 




86 88 



Figure 6: Numerical results for the method of Sec. 4. The two sides of the equation 
cos{BYA) = 27 as function of variable 7, including 30% SU(3) breaking effects in the 
p and a amplitudes. The physical value of 7 is determined by the intersection of the 
solid line with the wide band. The strong phases are taken as (0,0a) = (90°, 45°). 
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